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We present a model of neutrino oscillations in the framework of quantum field theory in which the 
propagating neutrino and the particles participating to the production and detection processes are 
described by wave packets. The neutrino state is a superposition of massive neutrino wave packets 
determined by the production process, as naturally expected from causality. We show that the 
(N . energies and momenta of the massive neutrino components relevant for neutrino oscillations are in 

' general different from the average energies and momenta of the propagating massive neutrino wave 

' packets, because of the effects of the detection process. Our results confirm the correctness of the 

standard expression for the oscillation length of extremely relativistic neutrinos and the existence 
' of a coherence length. 
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I. INTRODUCTION 



' The theory of neutrino oscillations is an important field of research, especially after the convincing evidence 

. of the existence of neutrino oscillations obtained in atmospheric and solar neutrino experiments (see [El 0, 

Neutrinos can oscillate if neutrino mixing is realized in nature (see ||^, ^ |10|, |11[), i.e. if the left-handed 
components i/ql of the flavor neutrino fields (a = e, ji, r) are superpositions of the left-handed components 
' VaL of the fields of neutrinos with definite mass: 



Qh! = 2^ Uaa VaL , (1-1) 



where U is the neutrino mixing matrix. Here the field Va describes a neutrino with mass rUa. A flavor neutrino 
- ^ , produced in a charged-current weak interaction process is a superposition of massive neutrino components 

■ with well defined kinematical properties. Neutrino oscillations are space-dependent and/or time-dependent 

Jh \ fiavor transitions generated by the different phase velocities of different massive neutrino components. The 

probabilities of flavor transitions in neutrino oscillation experiments depend on the differences Am^j, = 
— of the squared neutrino masses and on the elements of the mixing matrix U . 
Different models of neutrino oscillations have been presented by several authors (see [hO|, Jl2|, JlJ 
references therein). In this paper we adopt a wave packet approach, in the line of Refs. p5[ , Hfl, 

As discussed in Ref. , neutrino oscillations are possible only if the processes of neutrino production and 
detection have momentum uncertainties that allow the coherent production and detection of different massive 
neutrino components. From the fundamental principles of quantum mechanics it follows that the production 
and detection processes must be localized and the massive neutrino components of a flavor neutrino must be 
described by wave packets. We have developed a quantum-mechanical wave-packet description of neutrino 
oscillations in Refs. jl5[ This model provides important insights in the physics of neutrino oscillations. 
The main one is the confirmation of the correctness of the standard expression for the oscillation phase of 
extremely relativistic neutrinos (see 20, 2^, |2^). In addition, the quantum-mechanical model of neutrino 




oscillations implies the existence of a coherence length, which was predicted in Ref. pa on the basis of 
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intuitive quantum-mechanical arguments. For source-detector distances larger than the coherence length 
the flavor transition probability does not oscillate, because the separation of the different massive neutrino 
wave packets, which propagate with different group velocities, is so large that they cannot be absorbed 
coherently in the detection process |Q . In this case the probability of flavor-changing transitions is constant 
and depends only on the elements of the mixing matrix U. 

In the quantum-mechanical model of neutrino oscillations the expression of the state describing a flavor 
neutrino as a superposition of massive neutrino components has to be assumed, because quantum mechanics 
is not sufflcient for the description of the neutrino production and detection processes. Therefore, it is 
necessary to develop a description of neutrino oscillations in the framework of quantum field theory, which 
allows to calculate the amplitudes of neutrino production and detection. Moreover, one must notice that the 
theoretical prediction of neutrino oscillations follows from the mixing of neutrino fields in Eq. ( |l.l| ) and a 
consistent theory of neutrino mixing and oscillations can be formulated only in the framework of quantum 
field theory. 

In Refs. |l^ we have developed a quantum-field-theoretical model of neutrino oscillations in which 
the particles taking part to the neutrino production and detection processes are described by localized wave 
packets and the neutrino propagating between the production and detection processes is a virtual particle 

P8|, |29|| for another quantum-field-theoretical model 



(see also Refs. ||25|, 26 for similar approaches, Ref. 



of neutrino oscillations with virtual neutrinos, and Refs. |30, Q for a different quantum-field-theoretical 
point of view). This model confirms the correctness of the standard expression for the oscillation phase of 
extremely relativistic neutrinos and the existence of a coherence length. 

However, since in oscillation experiments neutrinos propagate over macroscopic or even astronomical 
distances, we think that it is unnatural to consider them as virtual particles, with undefined properties (see 
Refs. [HI 0). Since massive neutrinos propagate as free particles between production and detection, it 
should be possible to describe their superposition constituting the flavor neutrino created in the production 
process by an appropriate quantum-field-thcorctical state, similar to the quantum-mechanical wave-packet 
state in Ref. Causality demands that the neutrino state is determined by the production process. 

In this paper we present a quantum-field-theoretical model of neutrino oscillations in which the neutrino 
propagating between the production and detection processes is described by a wave packet state determined 
by the production process. This state is derived in Section starting from a production process in which 



the other interacting particles are described by wave packets. In Section [II we calculate the amplitude of the 



detection process occurring at a space-time distance {L,T) from the production process, using the neutrino 



state obtained in Section O and wave packets for the other interacting particles. In Section IV we calculate 
the probability of neutrino oscillations in space from the average over the unmeasured propagation time 
T of the detection probability and we discuss the effects of the detection process on neutrino oscillations. 
Conclusions are presented in Section 

In the following we assume that the particles that participate to the production and detection processes 
can be described by appropriate wave packets. This is possible if their properties are determined. If the 
information about their properties is incomplete, as often occurs in practice, each particle must be described 
by a statistical operator, also known as "density matrix" , constructed from an incoherent mixture of wave 
packets with definite properties. As a consequence, also the propagating neutrino must be described by a 
statistical operator constructed from an incoherent mixture of the pure states that we derive in Section ||, 
and the oscillation probability is given by an appropriate average of the oscillation probability that we derive 



in Section IV over the unknown properties of the particles participating to the production and detection 
processes. 

II. PRODUCTION 

The approach presented here is based on the fact that in quantum field theory the effect of interactions is 
described by the operator S — 1, with 



5 = T exp (^-i J d^xHiix)^ , 



(2.1) 
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where Tiiix) is the interaction Hamiltonian expressed in terms of the appropriate field operators. Given an 
asymptotic initial state \i) the asymptotic final state resulting from an interaction is given by 

I/) (X (5 - ^-i j d^xUiix) \i) , (2.2) 

where the last expression applies to first order in perturbation theory, which we adopt in the following, 
because we consider weak interaction processes. Hence, in quantum field theory it is possible to calculate 



with Eq.(2.2) the final state resulting from any interaction and in particular the final state of the processes in 
which a neutrino is produced. Since in the majority of experiments neutrinos are produced in charged-current 
weak decays, we consider the production process 

P/ ^ Pf + C + i^a , (2.3) 

in which P/ is the decaying particle, Pp is a decay product (absent in two-body decays) and £j is the final 
state charged lepton that determines the flavor of the produced neutrino Va- As a simple example of a 



possible production process of the type (2.2) we will consider the pion decay 

7r+ ^ /i+ + , (2.4) 

where P/ = 7r+ , Pp is absent and = /i+ 



Let us consider the production process (2.2). The final state obtained with Eq.(Ej) describes all the final 



particles of the process in an entangled way: 

|Pf,£+ I d''xn'i{x)\Pi), (2.5) 



where |P/) is the state describing the initial particle P/ and Ti^{x) is the Hamiltonian describing the 
production process. 

The entangled state \Pf, ^tn ^a) can be disentangled by measuring the properties of the particles involved. 
In particular, in the study of neutrino oscillations, one is i nterested in the knowledge of the state \ua) 



describing the neutrino produced in a process of type (2.3). In order to disentangle the state \i'a), it 
is necessary to measure the properties of the other particle in the final state, i.e. of Pp and t^. This 
measurement does not have to be done necessarily by a specific instrument, but could be done by the 
interactions of the anti-lepton £^ and of Pp with the surrounding medium. The measurement process causes 

a collapse of the entangled final state \Pp,t'^^Va) to a disentangled state \Pp)\i'^)\va) , with the neutrino 
state \va) given by 

\v^)^{{Pp\{l+J)\Pp,It,Vo.)^{{PFm\)-i fd^xnfix)\Pi). (2.6) 



The effective interaction Hamiltonian that describes the process (2.3) in the Standard Model is 



nfix) ^ ^p^{:,)jP{l-^,)Ux)J^{x) 

= ^Y.U*,D4x)r(.l-l5)iaix)J^{x), (2.7) 

where Gp is the Fermi constant and Jp{x) is the weak charged current that describes t he t ransition P/ Pp 



In order to simplify as much as possible our discussion, let us consider the process (2^) with the particles 
P/, Pp and described by the wave-packet states 

Ix) = j d'pi;^ip;p^,<JpMP^K))^ (2-8) 

where x — PiiPfi^^^ the momentum distributions are denoted by ip^{p;p^, ap^), and h-^ are the helicities. 
We assume that the particles P/, Pp and l'^ are polarized. If the measurement process is not sufficient 
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to determine the polarization of these particles, each of them must be descri bed by a statistical operator 
(density matrix) constructed from an incoherent mixture of the pure states (2^) with different helicities. 
Consequently, the propagating neutri no m ust be described by a statistical operator constructed from an 
incoherent mixture of the pure states {2£) obtained with different helicities of the particles Pi, Pp and 
We consider the Gaussian momentum distributions^ 



-3/4 



exp 



40-2 



(2.9) 



where p-^ and are, respectively, the average momentum and the momentum uncertainty of the wave 
packet of the particle x- The corresponding wave functions in coordinate space are given by 



ipx{x,t;p cTpx) = (0|x(a;)|x) 



(27r)3/ 



(2.10) 



where 



Exip) = \P 



(2.11) 



is the energy corresponding to the momentum p, a nd w e have neglected for simplicity the spin degrees 
of freedom. The Gaussian momentum distributions ( ^.9[ ) are assumed to be sharply peaked around their 
average momentum, i.e. the condition ^ E^{p-^)/mx is assumed to be satisfied. Under this condition, 
the energy E^{p) can be approximated by^ 



^xiP) - ^x+VxiP^Py) > 



(2.12) 



where 



ExiP 



x\i'x' 



(2.13) 



is the average energy (up to corrections of the order a'l /E^, see Eq. (2.17)) and 



dE-, 



dp 



(2.14) 



is th e gro up velocity of the wave packet of the particle x- Under this approximation the integration in 
Eq. (2.10) is Gaussian and leads to 



27ro-: 



-3/4 



xxl 



exp 



'iE^t + ip ■ X 



4a2 

"xx 



(2.15) 



where a^x defined by the relation 



'XX '-'PX 



(2.16) 



is the spatial width of the wave packet. One can see that the wave functions overlap at the origin of space- 
time coordinates, where the neutrino production process takes place. Hence, wave packet states (E^ are 
appropriate for the description of the particles taking part to the localized production process (2.3). 



^ A Gaussian momentum distribution is the most convenient one for the calculation of several integrations in the follow- 
ing. Other distributions which are sharply peaked around an average momentum p-^ lead to the same results after their 
approximation with a Gaussian in order to perform the integrals with a saddle-point approximation. Therefore, the Gaussian 
momentum distributions can be taken as approximations of the real momentum distributions from the beginning. 

^ With this approximation we neglect the spreading of the wave packets. 
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Let us determine the energy uncertainty of the wave packet state (2.8) describing a locahzed particle x- 
Developing E^(j}) up to second order in the power series of (p — p ), one obtains the average energy"^ 



{E)^ = {x\E\x)^E^+(3-v 



-2 \ 



2E^ 



where E is the energy operator. The average squared energy is given exactly by 

{E^)^ = {x\p'\x)+nil = El + 3al^, 

where P is the momentum operator, leading to the squared energy uncertainty 

{iSEr)^^{x\iE^E^f\x)^vlal^. 



(2.17) 



(2.18) 



(2.19) 



Therefore, somewhat surprisingly, a localized particle at rest has momentum uncertainty without energy 
uncertainty (at order ap^/E^). 

We use the following Fourier expansion for the spin 1/2 fermion fields: 



h=±i 



p'>=Ef{p) 



(2.20) 



where h is the helicity, Uj{p,h) and Vj(j)^h) are four-component spinors, and af{p,h) and bf{p,h) are the 
particle and anti-particle destruction operators obeying the canonical anticommutation relations 



{af (p, h),aUp', h')} - {b0, h),bUp', h')} - 5^{p - p') Shh' 



(2.21) 



The one-particle fermion states with definite momentum and helicity |f(p, h)) = a|(p, h) |0) are normalized 

by the relation (f(p, /i)|f(p', ft,')) — S^{p ~ p) 5hh' ■ One can easily check that in this way the wave-packet 

states (2.8) for y = fare nor mal ized to (f|f) — 1. 

From Eqs.(2^), (2/7) and (2^), the state that describes the neutrino produced in the process ( |2.3| ) is given 

by 

\va) cx ^u*^ d^p'p^ (ppp ; Ppj. , <7pPF ) / <i^p'i+ {Pt+ ; Pi+ . ) 

a J " " 

X J <i^PPj i^Pi {PPi ':PPj,CrpPi) 

X J d^x {{PF{p'p,)m{P,t)\) ^-i^)^' (1 - 75)^a(a;) j;(x) \Pi{p'p^)) . (2.22) 
Using the Fourier expansion ( ^.20| ) for the Icpton fields, we obtain 
Wa) ocJ^Ka J J d^p'p^^/j*p^{pp^;pp^,app^) 



(2.23) 



X J d^p[+ {p[+ ; p^+ , cr^+ ) J d^p'p^ ipp^ {pp^ ■,pp^, app^ ) {pp^ , hp^ ; pp^ , hp^ ) 
X / d'pY^u^Mhh'' il-l,)v,^{p[,,h,,Je'^'^''^t'^^^^^ 



I would like to thank M. Rpiithp fnr p oi ptintr o iit th p ner.es sity t n Hpv plop E-^{p) up to second order in the power series of 
(p — p^). As a consequence, Eqs. (2.17), ( |2.19| ), (2.52), (2.5S) and (2.70) have been corrected with respect to the fi rst v ersion 
of the paper appeared in the electronic archive hep-ph, where i?x(p) developed only to first order (as in Eq. ([2.12[)). 
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with / = E,,^{p), p'^^ = Epj{pp^), p'^^ = Eppip'p^), p'°+ = Eg+{p[+), and ((pp^,hpp;ppj,hp^) = 
{PF{ppp,hp^)\ Jp{0) \Pi{pp^,hpj)). For example, in the pion decay process ( |2^ ) Jp{pp^,hpj,;pp^,hp^) is 
given by (0| ^^(0) |7r+(p^+)) = /^r (p^+)p, where /tt is the pion decay constant. 

Since the wave packets of Pj , Pp and are assumed to be sharply peaked around the respective average 
momenta, the integrations over d^p'p^ , 'd^p'f + ^^id d^Pp^ can be performed with a saddle-point approximation 
using the approximation ( ^.12| ), which leads to 



J2uL [ d'pj2A^ip,h)w,{p,h)) 



X exp 



-i{Ep - Ey^ {p))t + i{pp - p)x - 



X -2vp ■xt + Y.pt^ 



with 



Ep — Epj — Epp — E^ 

Pp = Vpj - Ppp - Pi+ 
1 1 



^xPi 



VP = alp 



2 

^xP 



1 

1 


1 

1 




'2 2 
^xPp 


+ 


VP. 


«Pf , 






2 ~^ 
<^xP^ 




' -2 


-.2 
"Pf 


->2 
f ^ 


<P: 


^Ip^ 





-^a{P: h) = u^AP: h)Y (1 - 75) Jp{Ppp,hpp;p'p^,hpj) . 



(2.24) 

(2.25) 
(2.26) 

(2.27) 

(2.28) 

(2.29) 
(2.30) 



As naturally expected, the overall spatial width (TxP of the production process is dominated by the smallest 
among the spatial widths of P/, Pp and £'^. The quantities \vp\ and Sp are limited by 



0< \vp\ < 1, 0< Sp < 1. 
Carrying out the Gaussian integral over d^x, we obtain the neutrino state 

-sr(P) 



) = E / d'pe-'^^-^^ E -^aiP, h) h)) 



where Na is a normalization constant such that 



and 



Here 



S^aiP) 



{Va\Va) = 1 



_ (Pp - pf , [Ep - Ey^ (p) - {pp - p) ■ VpY 



4a^P 



+ 



A p = S p — w 



p 7 



such that 



< Ap < 1 , 

and we have defined the momentum uncertainty app through the relation 



axp cTpP — - . 



(2.31) 
(2.32) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 

(2.37) 
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The overall momentum uncertainty is dominated by the largest momentum uncertainty among P/, Pp and 



(2.38) 



The neutrino state (2.32) describes a neutrino produced in the weak interaction process (2.3) as a superpo- 
sition of massive neutrino components. In Section II] we discuss the detection of this state and in Section 
we derive the corresponding transition probability. In the following part of this Section we discuss some 
properties of the massive neutrino wave-packet states 



Wa)^NaJ d3pe-^r(?')^^r(p,/j)|j.,(p,;i)), 



(2.39) 



with the normalization constant Na such that 

{l^aWa)=l- (2.40) 

The massive neutrino wave-packet states ( ^.39 ) are the components of the state ( 2.32 ), which can be written 



as 



U* 



Na 

Li 

and the normalization constant Na is related to the normalization constants iV^ by 

-1/2 



(2.41) 



N^ 



UaaV 



1 X! 

\ a °- 

The normalization condition ( 2.40] ) requires that 



NlY, [ d'p\A^ip,h)\ 



2g-2Sr(p) 



(2.42) 



(2.43) 



The integration over d'^p can be done with a saddle point approximation around the stationary point of 



dp 



= 0. 



P=Pa 



The momentum p^ is given by 



Pp-Pa + ^ [Ep - Ea- {pp - Pa) ' V p] [Va - Vp) = , 



with 



EuAPa) 


\Jpl + m 


dE,M 


- P±. 


dp 





(2.44) 

(2.45) 

(2.46) 
(2.47) 



The saddle-point approximation of the integration over d'^p in Eq. (2.43) leads to 



N,. = 



flMAa 

\ 7r3 



1/4 



(2.48) 
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with 



A-' 



dpi dp'' 



, {^P <) (^^ - [{Ep - Ea) - (pp - p,) ■ vp] 6^'^ - vivt 



Hp 



2alp\p 



2alp\p 



E„. 



(2.49) 



Using the same saddle-point approximation for the integration over d^p, one can find that p^ is the average 
momentum of the state ji'a), 



{p)a = {Va\P\Va) = Pa , 

and the uncertainties of the three momentum components are given by 

((<5/)^)a = {'^aliP" - p'arWa) = \ {A-'f ■ 

Developing Ei,^ (p) up to second order in the power series of {p ~ p^) we obtain the average energy 



{E)a - {Va\E\Ua) ^ Ea + 

and the energy uncertainty 



^Ea 



(2.50) 



(2.51) 



(2.52) 



(2.53) 



In order to get furth er ins ight in the properties of the massive neutrino wave-packet states (2.39), it is 
necessary to solve Eq. ( 2.45| ) and determine the values of and Ea- It is important to notice that the 
production of the state \i/a) in Eq. ( ^.32 ) is not suppressed only if 

S^iPa) < 1 (2.54) 

for all values of the index a. Together with Eq. ( ^.45| ), this inequality constraints the possible values of p^, 
Pp and Ep. 

Equation (^.45 ) implies that the massive neutrino momenta must be aligned in the direction of pp: 



(2.55) 



with Pp = Pp £ a nd \£ \ = 1. Hence, £ is the unit vector in the direction of propagation of the neutrino. 

We solve Eq. ( 2.45| ) in the approximation of extremely relativistic neutrinos. This approximation is valid 
in practice because only neutrinos with energy larger than some fraction of MeV are detectable. Indeed, 
neutrinos are detected in: 

1. Charged-current or neutral-current weak processes which have an energy threshold larger than some 
fraction of MeV. This is due to the fact that in a scattering process 



(2.56) 



with A at rest, the squared center-of-mass energy s = 2E^mA + rn\ (neglecting the neutrino mass) 
must be bigger than {J2x leading to 



2mA 



rriA 
2 



(2.57) 



For example: 



9 



• El'^ ~ 0.233 MeV for i^e + ^^Ga ^ ^^00 + e" in the GALLEX ||, SAGE ||] and GNO [|| solar 
neutrino experiments. 

• E^J^ ~ 0.81 MeV for + ^''Cl ^'^Ai + in the Homestakc Q solar neutrino experiment. 

• E^ ~ 1.8 MeV for + p ^ n + in reactor neutrino experiments (for example Bugey p6|] and 
CHOOZ Hi). 

• El!^ ~ 2.2 MeV in the neutral-current process v + d^p + n + i' used in the SNO experiment to 
detect active solar neutrinos 

• El^ ~ llOMeV for + n ^ p + /.i- . 

• El^ ~ ml/2me ~ 10.9 GeV for i/^ + er -^Ve + M"- 



2. The elastic scattering process z^+e^ v+e^ , whose cross section is proportional to the neutrino energy 
{(j{Eu) ^ aQEy/rrie, with do ^ 10^^** cm^). Therefore, an energy threshold of some MeV's is needed 
in order to have a signal above the background. For example, E^^}^ ~ 5 MeV in the Super-Kamiokande 
[ p8[ solar neutrino experiment. 

Although the direct experimental upper limits for the effective neutrino masses in lepton decays are not 
very stringent {m^^ < 3eV, nii,^ < 190 keV, nii,^ < 18.2 MeV, see Ref. [^), we know that the sum of the 
masses of the neutrinos that have a substantial mixing with h'e , and v^- is constrained to be smaller than 
a few eV by their contribution to the total energy density of the Universe 0, |4^, . 

The comparison of the cosmological limit on neutrino masses with the energy threshold in the processes 
of neutrino detection implies that the main massive neutrino components of detectable flavor neutrinos are 
extremely relativistic^. 

We write the average massive neutrino energies Ea in the relativistic approximation'^ as 

Ea^E + C^, (2.58) 
where E is the neutrino energy in the limit of zero mass. The corresponding momentum has modulus 

p,^E-{1^0^- (2.59) 

In the following we consider 

Pp = Epi, (2.60) 

which corresponds to exact energy and momentum conservation in the production process in the case of 
massless neutrinos. It is possible to consider deviations from Eq. ( |2.60 ) compatible with the inequality (2.54), 



but such deviations would entail a considerable complication of the formalism without further insig 
the physical properties of neutrinos emitted in the process (|2.3|) 



its m 



* It is still possible that the three active flavor neutrinos t'e, t'/x, i't- have very small mixing with heavy massive neutrinos 
that are not relativistic in some experiment. In this case, the heavy neutrino masses must be taken into account in the 
calculation of the production and detection rates of the heavy massive neutrinos, but the oscillations generated by the large 
mass differences between light and heavy neutrinos are too fast to be observable. Therefore, in practice it is sufficient to 
consider an incoherent mixture of light and heavy massive neutrinos that generate a constant flavor-changing transition 
probability. Also the mass differences between possible heavy neutrinos are expected to be too large to generate observable 
oscillations. Therefore, in the following we study the oscillations due to the mixing of the flavor neutrinos i/^, and with 
light extremely relativistic massive neutrinos. 

^ In principle it is possible to consider the approximation of almost degenerate but not relativistic neutrinos . This 
approximation could be relevant in practice only if the three active flavor neutrinos Ve , Vfi, and Vt mix with heavy and almost 
degenerate massive neutrinos, such that the small mass differences among heavy neutrinos generate observable oscillations. 
Since this possibility seems very unlikely, we do not consider the case of almost degenerate but not relativistic massive 
neutrinos. Notice, howeve r t hat a similar approximation is important in the analogous quantum-field-theoretical treatment 
of meson oscillations (see [|13|). 
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At zeroth order in m'l/E'^, Eq. ( 2.45| ) imply 



E = Ef 



(2.61) 



The solution of Eq. ( 2.45 ) at first order in m'^/E'^ gives the value of ^: 

Xp -1 ■ vp (l ~ t ■ vp 



A, 



1 



■ Vp 



(2.62) 



Hence, the value of the parameter ^ that determines the average energy and momentum of the produced 
neutrino state depe n ds on the characteristics of the production process through the quantities Xp and vp. 



Using Eqs. ( |2.55| ), ( 2.60 ) and the relativistic approximations ( p. 5^ ), ( 2.59 ), with ^ given by Eq. (|2.62D , we 
find 



Xp+ l-I- 



Vp 



(2.63) 



Hence, the condition ( 2.54 ) is satisfied if 



l<'^ppE\ Xp+ (l-l-v 



(2.64) 



for all values of a. This inequality shows that a finite momentum uncertainty (Jpp of the production process 
is necessary in order to produce coherently different extremely re lativ istic massive neutrino components. 

The momentum and energy uncertainties in Eqs. ( 2.51 ) and ( 2.53| ) can be estimated at zeroth order in 
rrr^/E'^ inverting the symmetric matrix 



From the determinant 



we obtain® 



A = 



^ Ap Xp A" I 



V 



Ap 

(iip-l)t>f 



Ap 



1 



Ap ^ 

_ {vpY 



DctA = 



(2a^p)3 



(2.65) 



(2.66) 



{{5EY)a^-{K-') 



\\xx 



2 

<T.pP 



2 



1 + 



(^"p)'' + (^p)'' 
Ap 



Ap 
Ap 



2 



1 



Ap 

iv^p-lf + {vlf 



(t,g,_l)2 + (^^)2+(^. )2 : 



2 Ap 

^pp ^ {vf.-iy+{vir+{vpr ■ 



(2.67) 
(2.68) 
(2.69) 
(2.70) 



I would like to thank M. Beuthc for finding misprints in Eqs. (2.68) and (2.69) in the first version of the paper appeared in 
the electronic archive hep-ph. 
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Thus, the momentum and energy uncertainties are of th e order of cr^p, the total momentum uncertainty in 
the production process, unless the factors in Eqs. ( 2.67 )-( ^.70 ) depending on Xp and vp assume extreme 
values. 

In order to illustrate the meaning of our results, in the following Subsections we consider as an example 
r ieutr ino production in the pion decay process (2.4) at rest (1)7^+ = 0) in four cases. From Eqs. ( 2.6C| ) and 
( 2.61 ), the energy and momentum of the neutrino and muon at zeroth order in m^/i?^ are determined by 
energy-momentum conservation in the case of massless neutrinos: 



m I — m 



2m, 



29.8 MeV, 



2m, 



109.8 MeV, 



(2.71) 



leading to the muon velocity 



0.27. 



(2.72) 



Equation (2.19) implies that the localized pion at rest has no energy uncertainty (at order (Tp^+/m^+). 



A. Unlocalized production process 



In the limit 



0, a„ 







(7pP —^ , 



(2.73) 



the particles taking part to the production process and the production process itself are not localized. 

In this case the condition ( p. 64 ) is not satisfied for the coherent production of differen t ext remely rela- 
tivistic massive neut rino components is not valid. Furthermore, no deviat ion f rom Eq. (2.60) can satisfy 
the inequality ( 2.54 ) for more than one value of a, because in the limit ( 2.73 ) S^(p^) becomes infinite, 
suppressing the production of unless 



{Pp-Paf + ^[Ep~Ea 
Ap 



(Pp-Pa) -^P] 



0. 



(2.74) 



Taking into account the fact that Ap is positive, Eq. (2.74) can be satisfied only if both squares are zero, 
i.e. if pp — and Ep = Ea exactly. Since this constraint can be satisfied only for one value of the index 
a, only the corresponding massive neutrino is produced. 



B. Unlocalized pion 

If the momentum of the pion is determined with high accuracy, 

o'pT+ , 

the pion is unlocalized, 
and we have 

which imply (taking into account — 0) 



(2.75) 
(2.76) 
(2.77) 
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From Eq. (2.62) we get 



0.21 



(2.79) 



This is the value of ^ given by exact energy-momentum conservation in the production process in the case 
of pion decay at r est a nd different muon momenta for each massive neutrino. 
The condition ( 2.64 ) in this case reads 



< 



(2.80) 



If the muon propagates in normal matter, one can estimate (T^„+ to be given approximately by t he inter- 



10 eV. Thus, the condition ( 2.8C ) numeri- 



atomic distance, crx/i+ ~ 10 cm, which corresponds to (Tp^+ 

cally gives < 10^"'^ eV^, that is certainly satisfied. 

The production process has a mom entum unc ertain ty given by Eq. (|j7|), and an energy uncertainty 
Vfi+cFpP due to the muon. From Eqs. (|]67) and ( 2.70 ) one can see that in this case the massive neutrino 
components of the neutrino state ( ^.32| ) have no energy uncertainty and no momentum uncertainty along the 
direction £ of propagation. Hence, in this limiting case the energy and momentum uncertainties of the massive 
neutrino components of the neutrino state are rather different from the energy and momentum uncertainties 
of the production process. Only the unc ertai nties o f the components of the momentum orthogonal to the 
direction of propagation, given by Eqs. ( 2.68 ) and ( 2.69 ), are equal to the momentum uncertainty Upp of 
the production process. These uncertainties are actually necessary in order to localize the massive neutrino 
components along the direction of propagation. 



C. Equal energy limit 



In the limit in which the pion is localized but the final lepton is not localized. 



we have 



(2.81) 



<^xP 



which imply 



vp = Q.. 



(2.82) 



(2.83) 



In this limit the different massive neutrino components have the same energy. 



Since app = apT^+ ^ 0, the production process has a momentum uncertainty, but Eq. (2.1£) implies that 



both the pion and the muon have no ener gy un certainty. From Eqs. (2.67) and (2.70) it follows that each 
massive neutrino component in the state (2.32) has definite energy and momentum along the direction of 
prop agati on, without any spread. This is due to the fact that the relations ( 2.82|) imply that S^{p) in 
Eq. (2.34) is infinite, suppressing the production of i>a, unless E^^(j}) ~ Ep exactly. This constraint can 
be satisfied simultaneously for different values of the index a taking different values of p. The momentum 
uncertainty of the pion is necessary in order to allow the coherent production of a state with different massive 
neutrino components with different values of p. For each massive neutrino component there is only one energy 



Ea = Ep, equal for all components, and the corresponding momentum pa — yE' 



along the direction of propagation each massive neutrino component in the state (2.32) is a plane wave and 
not a wave p acket . 



This implies that 



Equations (2.68) and (2.6£) imply that the uncertainties of the components of the momentum orthogonal 
to the direction of propagation are equal to app. These uncertainties are allowed because they generate an 
energy uncertainty of higher order in Opp/Ep, that has been neglected in our formalism. Actually, as we 
have already remarked in the discussion of the previous example, these uncertainties are necessary in order 
to localize the massive neutrino components along the direction of propagation. 
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The case under consideration has some similarity with the one considered in Ref. [g7| by Grimus and 
Stockinger, in which it was assumed that the particles that take part in the production process (as well as 
those participating to the detection process) are in bound states with definite energy or are described by 
plane waves with definite energy. Although the physical picture in Ref. [2^ is different from the example 
under consideration, in which the localized pion at rest is free, in both cases there is no energy uncertainty 
in the production process and different massive neutrino components have the same energy because of exact 
energy conservation. In this sense, the case considered in Ref. pTf can be considered effectively as a limiting 
case of the general wave-packet treatment considered here, as done in Ref. |^ comparing the model of 
Ref. with the wave-packet model with virtual intermediate neutrinos discussed in Refs. ||l^, p^ . 

We think that in some cases it may be important to consider the fact that some particles are in bound 
states, but an appropriate description of such a case should take into account also the fact that the bound 
states (as atomic nuclei) are localized, leading to a wave packet description. We also think that the description 
of some particles taking part to the production (or detection) process by plane waves is unrealistic, because 
neutrinos are usually produced in dense media, where these particles are localized by interactions. 

In conclusion of this subsection we would like to remark that, although the equal energy limit that we 
have considered is realizable in principle, we think that in practice it is rather unlikely since the pion (or in 
general the initial particle Pi) must decay exactly at rest and the produced muon (or in general the final 
particles Pp and £j) must be completely unlocalized. If the pion does not decay at rest in the reference 
frame of the observer, it is possible to boost the reference frame to the one in which the pion is at rest, but 
the pion velocity in the reference frame of the observer must be known with high accuracy. This information 
is usually not available, for example in the existing accelerator and atmospheric neutrino experiments in 
which neutrinos are produced by pion decay in flight. 



D. Realistic case 



In a realistic experimental setup the localizations of the pion and muon are of the same order of magnitude. 
Indeed, the typical neutrino production in pion decay at rest^ occurs in a medium, where both the pion and 
muon are localized by interactions with the surrounding atoms. Let us consider, for example. 



which leads to 



0'a:7r+ — f^xp+ — 2CTj;P , (2.84) 



wp ~ — - —I, i.p~-— , Ap~— — . (2.85) 



2 2 ' 2 ' 4 

In this case we have 



1 + 



0.13. (2.86) 



Hence we see that the value of ^ in a realistic situation is of the same order of magnitude as that in Eq. ( 2.79 ), 
corresponding to exact energy-momentum conservation in the production process in the case of pion decay 
at rest and different muon momenta for each massive neutrino. 



In accelerator experiments in which neutrinos are produced by pion decay in flight the localization of the pion and muon 
are given by the dimensions of the decay tunnel. In solar neutrino experiments electron neutrinos are produced in the core of 
the sun where all the particles taking part to the production process are localized by interactions with the dense medium. In 
atmospheric neutrino experiments neutrinos are produced by pion and muon decay in flight in the atmosphere, where pions, 
muons and electrons are localized by the interactions with air. In reactor neutrinos experiments electron antineutrinos are 
produced by the decays of heavy elements in the dense reactor core where the heavy nuclei and electrons are localized by 
interactions with the medium. 
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The condition (2.64) in this case becomes 



< 



(TpP E ' 



2 



E. 



(2.87) 



If the pion decay occurs in normal matter, UxP is given ap proxi mately by the inter-atomic distance, <TxP 
10~^ cm, corresponding to Opp 
certainly satisfied. 



10^ eV, the condition ( 2.87 ) numerically reads < 10^ eV^, that is 



Since in this realistic case vp and Ap do not have extreme values, from Eqs. ( 2.67 )-( 2.70 ) one can see that 
the energy and momentum uncertainties of the neutrino state are of the order of <Jpp, about 10^ eV for pion 
decay in normal matter. 



III. DETECTION 



Let us consider neutrino detection through the charged-current weak process 

vp + Di^Dp+l-p, (3.1) 

at a space-time distance { L, T) from the production process. 

The state \va) in Eq. (2.32) describes the neutrino produced in the process (2.3) at the origin of the 
space-time coordinates. Since we want to describe with the same formalism the detection of the neutrino 
through the process (3.1) occurring at a space-time distance {L,T) from the production process, we must 
translate the origin to the detection space-time point. Hence, the neutrino state relevant for the detection 



process is obtained by acting on \ va) with the space-time translation operator exp y—iET + iP • L j , where 
E and P are the energy and momentum operators, respectively. The resulting state is 

\v^{L, T)) =N^J2 / exp (-lE,^ {p)T + tp ■ l) e-^"(?) ^ A^{p, h) \va{p, h)) . 



(3.2) 



The amplitude of interaction of the neutrino state \i'a{L^T)) in the detection process (3.1) is 
Ao,p{L,T) = {DF,rp\ - i Jd^xnfix) \Di,y^{L,T)) , 

where 



(3.3) 



Gp 
G 



^ Upt ip{x) Y (1 - 75) ^b{x) J^(x) , 



(3.4) 



is the effective interaction Hamiltonian that describes the detection process, and J^{x) is the weak charged 
current that describes the transition Di Dp. 

For simplicity, we assume that the particles that take part to the detection process are described by the 
Gaussian wave-packet states (2^) with x — Di^Dp^l'^. Using the neutrino state (3^) and the Fourier 
expansion ( 2.2Cl| ) for the lepton fields, respectively, we obtain 



I e ^^^^P^ exp 



A^p{tT)^Y.^lJJpa d^xY, / d?pAP,{p,h) 

a h 

X j d^p'c^ i^*Dp (Pdp ; Pdp , f^pDp ) J d^Pf- ^l- {p'l,- ; Pi,- , cr^- ) 

X / <i^P'D,'^Di{PDi'^PDi^(^pDi) 



iEy^ {p)T + ip- L 



nil iip'Dp+P'i,-~P'Dj-P)^ 
^u^^^{p^-,h^-)Y {l~l5)u^AP-,h)Jp {PDF,hDp]PD,-,hDi)e f , (3.5) 
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where / = E^^{p), where p'°^ = EdAp'di)^ p'dj. = Edj,{pd^), pf- = E^-{p[-), and 

Jpip'op^^Dp^PDi^hD,) = {DF{PDp,hDp)\ Jp{0) \Di{p'oj,hDi))- Following the same method as that used 
for the production process (see Eq. ( 2.24 )), we perform the integrals over d^p^^, d'^p^_ and d^p'jj^ with a 

saddle-point approximation through the approximation (2.12), obtaining 

Ao^pitT) (X Y.^lJ}paY. I d'pA^ip,h)A^ip,h)e-'^"^P^ exp 



■'iEi,^{p)T + ip ■ L 



X exp 



-i{E^^{p) - ED)t + i{p-pjy)x - 



2vD ■ xt + Y^ot^ 



with 



Er 



Er 



Ef- — Er 



Pd = PDj. + Pi- - Pd, 




= CTxD 



{P, h) = ug-{pg- , hf-) Y - 15) (p, h) {pjj^ , ho J, ; Pdj , ho J ) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 

(3.11) 
(3.12) 



As expected, the overall spatial width Uxd of the production process is dominated by the smallest among 
the spatial widths of the participating particles. As the corresponding quantities in the production process, 
\vd\ and are limited by 



0< < 1, 0< Ed < 1. 
The Gaussian integral over d'^x leads to 

Ao.p{L, T)^Y. U^aUpa I d'p-4r(p, h)A^{p, h) e-^"(^) exp 



-iE^Sp)T + ip ■ L 



where 



Sa{p)=S^{p)+S^{p) 



(3.13) 



(3.14) 



(3.15) 



and S^{p) has the same structure as S^{p) given in Eq. ( 2.34 ), with the quantities relative to the production 
process replaced by the corresponding ones relative to the detection process: 



S^ip) 



_ {Pd - Pf , [{Ed - E^^ (p)) - {p^ - p) ■ vdY 



4 



4^2^ Ad 



with 



limited by 



Ad — Sd — Vjj 



< Ad < 1 , 



(3.16) 

(3.17) 
(3.18) 
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and the momentum uncertainty defined by 



(3.19) 



which is dominated by the largest of the momentum uncertainties of Dj, Dp and 



2 2 



pDi 



'pDf 



(3.20) 



In general the dominant momentum contribution to the integration over d^p ( ^.14| ) does not come from 
the stationary point of S^{p), but from the stationary point of Sa{p), which takes into account also the 
momentum and energy uncertainties of the detection process. If these uncertainties are smaller than those 
of the production process, the detection process picks up as dominant co ntr ibution to the flavor-changing 
amplitude (3.14) a value of the neutrino momentum in the wave packet (3.2) that is significantly different 
from corresponding to the stationary point of {p) . 

We denote by ka the momentum corresponding to the stationary point of Sa (p) , 



dp 



= 0, 



(3.21) 



P—^a 



which gives the dominant contribution to the transition amphtude (3.14). The value of ka is given by 



Pp ~ka Ep^^a- [Pp -ka)-Vp 



'pP 



oIp Xp 



+ 



Pd -ka ^D-ea- [Pd - j ' «I? 



■'pD 



with 



Sa = E^^{ka) = V k„ 



[Ua - Vp) 



{Ua -Vd) =0, 



'a 1 



dp 



P—ka 



(3.22) 

(3.23) 
(3.24) 



We have chosen the notation fc^, and Ua in order to emphasize that in general these quantities are not 
the average momenta, energies and velocities of the neutrino wave packets propagating between production 
and detection, denoted by p^, Ea and Va, which are determined only by the production process, as follows 



from Eqs. ( 2.45 )-( 2.47 ). The quantities ka, Sa and Ua are approximately equal to the average neutrino wave 
packets momenta, energies and velocities p^, Ea and Va only if a^p, 3> Upp and a'^p, Xp ^ o'pp Xp. In other 
words, the properties of the neutrino wave packets can be measured only with a detection process having 
a relatively large momentum and energy uncertainty, which correspond to a relatively sharp spatial and 
temporal localization. 



Before solving Eq. (3.22), we notice that the amplitude (3.14) is not suppressed only if 



Saika) < 1 



(3.25) 



for all values of a. The inequality ( 3.25| ), together with Eq. ( 3.22 ), constraint the possible values of ka, 
Pp, Ep, Pp and Ep. Since both S'^(fca) and Sa{ka) are positive, we have the conditions S'^(fca) <C 1 and 
Saika) ^ 1 for all values of a. The first inequality is similar to the condition ( ^.54 ), but now it concerns 
the momenta ka and energies Sa that are relevant for the flavor transition amplitude (3.14). 

Similarly to what we have done in the discussion of the production process (see Eq. ( 2.60D ), we consider 



Pp ^ Ep' 



Pd^ Ept 



Ep — Er 



(3.26) 
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which corresponds to exact energy and momentum conservation in the production and detection processes in 
the case of massless neutrinos. Here £ = L/|L| is the unit vector in the direction of propagation of the neutrino 
from the production to the detection processes. It is possible to consider deviations from the relations 
(3.26) compatible with the inequality ( 3.25| ), but such deviations would introduce many complications in the 
following formalism, without further insights in the physics of neutrino oscillations. 

Equation ( |3.22 ) implies that the massive neutrino momenta ka must be aligned in the £ direction: 



ka — ka £ • 



We solve Eq. (|3.22|) in the relativistic approximation 



(3.27) 



At zeroth order in rri^/E'^ we obtain 



9 

771 

ka^E-{l~p)^ 



E = Ep = Ed , 



(3.28) 
(3.29) 

(3.30) 



and at first order 



l-vp[l~-l-vp) I-vdI^I-I-vd) 



p = 



>p - 



1 {1-1-vpY (i-vn-iy 



with 



1 



(3.31) 



(3.32) 



Notice t hat i n the pro duction and detection processes the squared momentum uncertainties add, as shown 
in Eqs. (2.27) and (3.9|), whereas in the total amplitude the inverses of the squared momentum uncertainties 
of the production and detection processes add. This is expected on the basis of simple physical arguments. 
Indeed, a large momentum uncertainty of a particle must increase the total momentum uncertainty in the 
corresponding process, whereas a small momentum uncertainty in one of the two processes constraints the 
momentum uncertainty of the neutrino connecting the two processes leading to a restriction of the momentum 
uncertainty in the other process. On the other hand, in th e pro duction and detection processes the inverses 
of the squared space uncertainties add (see Eqs. ( 2.27 ) and (O)), whereas in the total amplitude the squared 
space uncertainties of the production and detection processes add: 



(3.33) 



with ax defined by the relation 



'X 



(3.34) 



Also the behavior in Eq. ( 3.33| ) is expected on the basis of simple physical arguments: a small space uncer- 
tainty of a particle localizes better the corresponding process, whereas a large space uncertainty of one of 
the two processes increases the coherence of the overall process. 
In the relativistic approximation Sa{ka) is given by 



(- 



K'iEa, 



(3.35) 
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with 



c = 



i~e-vp 



[l-ivp-vo)] 



Ap cr-f Ar) 



(l~2-5p) (1-2-5d 



It follows that in the relativistic approximation the condition (3.25) reads 



(3.36) 



(3.37) 



If this inequality is satisfied for all values of the index a, the choices (3.26) are acceptable and the different 
extremely relativistic massive neutrino components contribute coherently to the flavor-changing transition 
amplitude. 



Let us return to the calculation of the transition amplitude. The integral over d'^p in Eq. ( 3.14 ) can be 
performed with a saddle-point approximation around the stationary point of (p), leading to 



■^ajka, h)Aa{ka, h) _Sa{ka) 



X exp 



„r + .fc, . i - i {L^ <T) {n-y^ {L^ - u^T) 



(3.38) 



where 



2^ 



p—k 

.3 



{Ep - Ea) - (pp - fca) 



Vp 



6^^ 



li Ilk 



2<Jlp\p 



2a2pAp 



{Ed - £a) 



Pd 



vd 



6^ 



The factors 



(3.39) 



(3.40) 



in the exponential of Eq. ( 3.38| ) do not suppress the transition amplitude only if L is aligned with the 
common direction i of the momenta ka, apart from a deviation of the order of fla/L, which is very small 
if the localizations of the production and detection processes are much smaller than the source-detector 
distance, a condition which is necessary for the observation of neutrino oscillations and which is satisfied in 
all neutrino oscillation experiments. Therefore, we consider L ~ LI and we align i along the direction of the 
X axis. Thus, the amplitude (3.38) can be written as 



A^,{L, T) oc ^ U:^Upa E e-^°(^°) exp 



VDetfia 



-i£aT + ikaL 



where 



(3.41) 



(3.42) 



are the spatial coherence widths. 
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In the relativistic approximation the product A^(ka,+)A^(ka,+) corresponding to the positive hehc- 
ity component of the massive neutrino Va is suppressed by the ratio m'^/E'^ with respect to the product 
A^ika, —)Aa{ka, — ) corresponding to the negative hehcity component and can be neglected. In the same 
approximation, the factors 



AP{ka,~)A^{ka,-) ^_s^Ck.) 



(3.43) 



can be approximated with their value in the case of massless neutrinos, can be extracted out of the sum 
over the index a and absorbed in the overall normalization factor of the flavor transition amplitude. The 
factorization of these quantities allows the calculation of the flavor transition amplitude independently from 
the production and detection rates. The oscillation probability obtained from this flavor transition amplitude 
can be used in the usual calculations of event rates in neutrino oscillation experiments, given by the product 
of the neutrino flux calculated for massless neutrinos, times the oscillation probability, times the detection 
cross section calculated for massless neutrinos (see |l^, [ll|]). 

It is important to notice, however, that a special care is needed for the factors e~^''^'''^\ to make sure 
that they do not suppress the contribution of the corresponding massive neutrino component. The factors 
g-Sa{ka) QQ^j^ ]-,g Qg^j^ ]-,g approximated with their value in the case of massless neutrinos only if 



Sa{ka) « 1 

for all values of the index a. Hence, the condition ( |3.37 ) must be strengthened to 



(3.44) 



(3.45) 



for all values of a. Let us emphasize that this condition is necessary for the unsuppressed production and 
detection of the different extremely relativistic massive neutrino components whose interference generates 
neutrino oscillations. In principle, it is possible to consider degenerate neutrinos for which the massless 
approximation is not appropriate, but as noted in footnote this case is irrelevant in practice. In any case, 
considering a case in which the massless approximation is not appropriate would require the inclusion of the 
neutrino mass effec t, wh ich are normally neglected, in the calculation of the production and detection rates. 

If the condition (3.4?;) is satisfied for all values of a, the transition amplitude in the relativistic approxi- 
mation can be written as 



-iSaT + ikaL 



[L - UaTf 



4^2 



Acf3{L, T)(xY^ U*^Uf3a exp 

a 

Here Sa and ka are given by their relativistic approximations ( ^3.28 ) and ( ^.29 ), and 

Ua ^ 1 



2E^ 



(3.46) 



(3.47) 



The coherence widths rja have been approximated with their value ry at zeroth order in ml/E"^. 

The value of ry is giv en by the inversion of the zeroth order approximation in ral/E"^ of the symmetric 
matrix Q.a in Eq. ( 3.39| ), 



2ai 



2alpXp 



O 



E^ 



(3.48) 



We obtained 



2 2 



(3.49) 
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with 



UJ 



l)vl 



{vl - 1) vlf + [{v-p - 1) t^f, - K - 1) v^, 



^Ip^^p 



'pD 



'pP 



Xp (Tpp,\ 



(3.50) 



Let us notice that the width rj obtained in Eq. ( ^.4£ ) is somewhat different from the corresponding width 
obtained in Ref. iQ. The differ ences stem from the fact that J7a in Eq. ( 3.3E ) is a matrix, because of the 
integration over d^ in Eq. (3.14), whereas the width in Eq. (18) of Ref. |18 is a number obtained from the 
integration over dq° in Eq. (15) of Ref. |p^ . 

In the hmit (ipp <^ apo the total momentum uncertaint y gp a nd the coherence width rj are dominated 
by the production process, and p ~ ^, with ^ given in Eq. ( 2.62| ). This happens if the production process 
is much less localized than the detection process. In this case ka — Pa, the average momentum of the 
massive neutrino component Va of the state {va) created in the production process. In other words, if the 
detection process is much more localized than the production process, i.e. if the momentum uncertai nty o f 
the detection process is much larger than that of the production process, the transition amplitude ( 3.46| ) 
depends only on the properties of the neutrino created in the production process. 

However, in general, Up, 77 and p depend on both the production and detection processes and in the limit 
o'pD o'pP, in which the detection process is much less localized than the production process, they are 
dominated by the detection process. 

The present derivation of the transition amplitude in neutrino oscillation experiments is more complete 
than the simple quantum-mechanical model presented in Refs. jl^, 16 , in which the energies and momenta 
of the massive neutrino components contributing to the transition amplitude are undetermined and have to 
be assumed. 



Formally, the transition amplitude (3.14) can be derived by projecting the state \va{L,T)) in Eq. (3.2) on 
the state 



a •' h 

representing the detection process: 

Aa0{L,T) CX (iyp\VaiL,T)) , 



(3.51) 



(3.52) 



as in the quantum-mechanical derivation of the flavor transition probability in neutrino oscillation experi- 
ments. However, the calculation of the coefficients of the massive neutrino components in the states {L, T)) 
and \1yf3) require the quantum field theoretical framework that we have adopted. 



IV. TRANSITION PROBABILITY 



In this Section we calculate and discuss the transition probability in space obtained from the average of 

P„0(i,T)cx |^„0(L,T)|2 (4.1) 

over the unmeasured propagation time T, as done in Refs. [|l^ |l^, |l^, Let us notice, however, that 
although in present neutrino oscillation experiments only the source-detector distance L is known, in the 
future it may be possible to measure also the propagation time T, leading to the experimental relevance of 



the space-time dependent transition probability (4.1) (see Refs. |43, 44 
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From the Gaussian integration over dT of Pai3{L,T), in the rclativistic approximation we obtain^ 

2 



a>b 



-2TTi- 



L 

T OSC 



T coll 

^ab 



T OSC 

^ab 



(4.2) 



where 



are the oscillation lengths, and 



T OSC 

^ab 



T coll 

^ab — 



\Ar 



(4.3) 



(4.4) 



are the coherence lengths. 

The transition probability (4.2 ) has the same form as that obtained in Ref. ||l^ in quantum field theory with 
virtual propagating neutrinos. The value of p in quantum field theo ry with virtual propagating neutrinos, 
derived in Ref. is the same as the one derived here (see Eq. (3.31)). Instead, the value of to is somewhat 
different from that derived in Ref. as already remarked after Eq. (3.5C). 

The last term in the exponential of the transition probability (4.2) suppresses the corresponding oscillatory 
term unless the localization of the production and detection processes is much smaller than the oscillation 
length, 



OSC 

b ■ 



(4.5) 



This is due to the average of the transition probability (4.1) over the unmeasured propagation time T. 
Indeed, a spatial uncertainty (j^ hiiply a similar time uncertainty, as one can understand from the relation 
between the energy and momentum uncertainties in Eq. ( 2.19| ), or by noticing that the spatial region in 
which a process can occur coherently must be causally connected. If the time uncertainty, of the order of 
ax , is larger than the oscillation length, the average of an oscillatory term over the propagation time T is an 
average over all oscillation phases which depend on T, a nd the result is zero. 

However, as discussed in Ref. the condition (4^) for the observability of neutrino oscillations is not 
necessary if 



0. 



(4.6) 



As shown in the example presented in Subsection IV B. this can be obtained with p = 0. In this case, 
since the different massive neutrino components have the same energy Sa = E, the oscillation phases do 
not depend on T and the average over T of the space-time dependen t tr ansition probability ([4T|) has no 



effect. Furthermore, in this case the coherence lengths LJj'^'^ in Eq. (4^) can be increased without limit 



| p6| . However, one must notice that this unlimited increase must be obtained by increasing uj and not , 
because an unlimited increase of ax would bring ap and a violation of the condition ( |3.45 ) necessary 
for the coherent production of the different massive neutrino components whose interference generates the 
oscillations. 



In the following Subsections we consider three cases analogous to those considered in Subsections [I A , 



lie and II D, that illustrate the effects of the detection process. In all these cases we consider the initial 



detection particle D/ at rest {vdi — 0). 



We order the massive neutrinos by increasing mass: mi < m2 < . . . . In this way the Am'^^'s in the sum over a > b are all 
positive. 



22 



A. Unlocalized detection process 



If the particles taking part to the detection process are unlocahzed, we have the hmit 



(4.7) 



In this case the condition ( p. 371) is no t satisfied. Moreover, si milar ly to the case of an unlocali zed p roduction 
process discussed in Subsection II A, no deviation from Eq. (3.26) can satisfy the inequality (3.25) for more 

than one value of a. Indeed, in the limit ( [f.7| ) S^ikg) becomes infinite and suppresses the production of Va 
unless 



[Pd 



ka\ + — 



1 



Ed 



(pD -ka'j 



VD 



0. 



(4.8) 



Since Xd is positive, Eq. (4.8) can be satisfied only if both squares are zero, i.e. if = ka and Ejj — 
exactly. Since this constraint can be satisfied only for one value of the index a, only the corresponding 
massive neutrino is detected and there are no oscillations, which are due to the interference of different 
massive neutrino components. 



B. Equal energy limit 



Let us consider a localized initial detection particle Dj at rest, and unlocalized final detection particles, 



0, 



which imply 



'pDi , 



= 0, 



In this case, Eq. (2.19) implies that all the detection particles have no energy uncertainty. 
Equation ( |3.31 ) gives 



p = 0. 



(4.9) 



(4.10) 



(4.11) 



Hence, the detection process picks up the momenta ka of the different massive neutrino components of 



the state \va{L,T)) in Eq. (|3.2D corresponding to the same energy Sa = E, corresponding to exact energy 
conservation in the detection process. These momenta are different from the average momenta of the 



massive neutrino wave packets that constitute the state \va{L,T)) , which are given by Eq. (2.4£), unless the 
conditions of Subsection [I C are satisfied by the production process. 
From Eq. (|3.36D we find 



4 {^-"^^y 



Thus, the condition ( p. 45 ) can be satisfied by a sufficiently large ap. 
Since the values in Eq. (4.1C) imply that 



in the case under consideration 



r coh 



+00 , 



+CXD , 



(4.12) 



(4.13) 



(4.14) 



i.e. the oscillations remain coherent at arbitrarily large distances. In agreement with the arguments presented 
in Ref. |2fl], the last term in the exponential of Eq. (4.2) does not suppress the oscillations because 



(4.15) 
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Physically, the infinity of the coherence lengths is due to the fact that the detection process picks up the 
momentum components with equal energy of the different massive neutrino components. In this case there 
arc no oscillations in time and the average over time of the space-time dependent transition probability is 
irrelevant. 

Notice, however, that the coherence lengths are i nfinit e because of the infinity of co, whereas the uncertainty 
ax has to be finite in order to satisfy the condition ( 3.45 ) necessary for the coherent production and detection 
of the different massive neutrino comp onen ts. 

As we have remarked in Subsection II C for the similar limit in the case of the production process, the 
required conditions {vdi — and Eq. (4.9)) are unlikely to be achieved in any realistic experiment. We think 
that in the case of the detection process the required conditions are even more unlikely to be achieved than 
in the case of the production process, because neutrino detection requires in practice to reveal at least one 
of the final pa rticl es in the detection process. This implies that these particles cannot be unlocalized and 
the condition (4.9) is unrealistic. 



C. Realistic case 



In a realistic case production and detection occur in matter and all the detection particles have uncertainties 
of the same order of magnitude, as well as all the production particles, as we have discussed in Subsection II D| . 
Let us consider, for example, the case in which the order of magnitude of the spatial localization of the 
detection particles is much larger than that of the production particles: 



O'xDi — CFxDf 



(4.16) 



which imply that the momentum uncertainty of the detection process is much smaller than that of the 
production process, and 



-.2 , -.2 



(4.17) 



Let us consider, for example, a muon neutrino (a = /i) produced in the pion decay process (2.4) at rest 
and detected as an electron neutrino (/3 = e) in the process 



g.s. 



(4.18) 



where ^^Ng.g. is the ground state of ^^N, with an atomic mass excess AM(^^N) ~ 17.3 MeV. Neglecting the 
nuclear recoil energy, the energy of the electron is given by 



E^- ~E- AM('^N) ~ 12.5 MeV, 



(4.19) 



where E is the energy of a massless neutrino given in Eq. ( |2.71 ) . Let us consider a recoil electron emitted 
in the backward direction: 



(4.20) 



Since \vi 



\'"Dt 



= and 



1, we have 



From Eq. ( 3.31 ) we get 



VD 



0.29. 



(4.21) 



(4.22) 



This value of p is different from the value of ^ in Eq. ( 2.8(: ), albeit of the same order of magnitude. There fore, 
the energies and momenta of the massive neutrino contributions to the ffavor transition amplitude (3.46) are 



24 



different from the average energies and momenta of the massive neutrino components of the neutrino state 
created in th e pro duction process. 
From Eq. (3.50) we obtain 



17, 



and from Eq. (4.4) we obtain the coherence lengths 



2 X IQis eV^ 



''ah 



■ CTxD ■ 



(4.23) 



(4.24) 



If jAm^jjl ^ 1 eV^ and a^n is given approximately by the inter-atomic distance, (Tj;jy ^ 10 ^ cm, we have 
i™'' ^ 10^ k m, a sufficiently long coherence distance for short-baseline neutrino oscillation experiments. 



From Eq. ( 3.36 ) we find that is approximately given by 

Xd 



c 



Xn+il-i-VD 



0.53. 



and the condition (3.45) is satisfied for 



(4.25) 



(4.26) 



If ffxD is given approximately by the inter-atomic distance, we have apo ~ 10 eV and the inequality (2.54) 
reads ^ lO-'^^eV^, that is certainly satisfied. 



V. CONCLUSIONS 



We have presented a quantum-field-theoretical model of neutrino oscillations in which the neutrino prop- 
agating between the production and detection processes is described by the wave-packet state (2.32), which 
is determined by the production process as naturally expected from causality. Since in real oscillation exper- 
iments neutrinos propagate over macroscopically large distances, we think that this model is preferable over 
the quantum-field-theoretical model of neutrino oscillations with virtual intermediate neutrinos presented in 
Refs. 00. 

We have considered production and detection processes of the form (|2.3|) and (3.1) (other types of pro- 



duction and detection processes can be considered with straightforward changes to the formalism) in which 
the interacting particles are described by localized wave packets. The Gaussian form of the wave packets of 
interacting particles that we have adopted should be considered as an approximation of the real wave pack- 
ets. In any case, other wave packets which are sharply peaked around their average momentum lead to the 
same results, because these results depend on several integrations that are calculated with the saddle-point 
approximation. 

We have obtained the flavor transition probability, presented in Eq. (4.2), that is almost identical to 
the one derived in the q uantu m-field-theoretical model with virtual intermediate neutrinos p^ . Only the 
quantity w, given in Eq. (3.50), is somewhat different from the corresponding one in Ref. JlSt . 

As already discussed in Refs. |2^, the flavor transition probability is determined not only by the 
production process, but also by the detection process [ p4[ . In particular, we have shown that the energies 
and momenta of the massive neutrino components relevant for the oscillations are in general different from 
the average energies and momenta of the massive neutrino components of the propagating neutrino state, 
which are determined only by the production process. 

Our result confirms the correctness of the standard expression ( |4.3| ) for th e os cillation lengths of extremely 
relativistic neutrinos and the existence of coherence lengths given by Eq. (44). We agree with the author 
of Ref. |2^ on the possibility to extend without limit the coherence length with an appropriate setup 
of the detection process. We have shown that the coherence length cannot be increased without limit 
by decreasing the momentum uncertainty of the detection (or production) process, because a vanishing 
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momentum uncertainty does not allow the coherent detection (or production) of different massive neutrino 
components. Instead, the coherence length ca n be increased without limit by cho osing a detection (or 
production) process such that p, given in Eq. ( 3.31 ), is zero and u, given in Eq. ( |3.50| ), is infinite. As 



discussed for the example presented in Subsection IV B , in practice the possibility to increase without limit 
the coherence length is rather unrealistic. As shown in Subsection IV C, in a realistic experimental setup 
the coherence length is very long, but not infinite. 

Finally, let us recall the important remark presented at the end of the introductory Section I. In th is p aper 
we have assumed that the particles Pi, Pp and participating to the neutrino production process (2.3) are 



described by the pure wave-packet states (2.8) in which all their properties are determined. In practice it is 
common that the knowledge of these properties is less than complete. In this case the particles P/, Pp and 
£j must be described by statistical operators (density matrices) constructed from appropriate incoherent 
mixtures of the pure wave-packet states (2.8). Consequently, the neutrino created in the production process 



must b e des cribed by a statistical operator constructed from an incoherent mixture of the pure wave-packet 
states ( ^.32 ). Similarly, if there is incomplete knowledge of the properties of the particles P/, Pp and 
and the particles Z?/, Dp and participating, respectively, to the production and detection processes, the 
oscillation probability is given by an appropriate average of the probability in Eq. (4.2) over the unknown 
quantities. 
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